Five Dimensional Cosmological Models in 
General Relativity 

R. K. Tiwari 1 . Farook Rahaman 2 . Saibal Ray 3 

o: 

Abstract 

o; 

A Five dimensional Kaluza-Klein space-time is considered in the presence of a per- 
fect fluid source with variable G and A. An expanding universe is found by using 
a relation between the metric potential and an equation of state. The gravitational 
constant is found to decrease with time as G ~ £~( 1-w ) whereas the variation for 
the cosmological constant follows as A ~ £~ 2 , A ~ (R/R) 2 and A ~ R/R where u 
1 fj 1 ' is the equation of state parameter and R is the scale factor. 
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1 Introduction 

o: 

>• . In recent years there has been considerable interest in the cosmological mod- 

els with variable gravitational constant G and the cosmological constant A. 
Variation of the gravitational constant was first suggested by Dirac [1] in an 
attempt to understand the appearance of certain very large numbers, when 
atomic and cosmic worlds are compared. He postulates that the gravitational 
constant G decreases inversely with cosmic time. Canuto et al. [2113] made 
numerous suggestions based on different arguments that G is indeed time de- 
pendent. Beesham [3] has studied the creation with variable G and pointed 
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out the variation of the form G ~ i -1 , originally proposed by Dirac [lj. On 
the other hand, Einstein introduced the cosmological constant A to account 
for a stable static universe as appeared to him at the time. When he later 
knew of the universal expansion he regretted its inclusion in his field equa- 
tions. Now cosmologist believe that is not identically, but very close to zero. 
They relate this constant to the vacuum energy that first inflated a universe 
causing it to expand [5] . From the point of view of particle physics a vacuum 
energy could correspond to quantum field that is diluted to its present small 
value. However, other cosmologists dictate a time variation of this constant 
in order to account for its present smallness [BJ. The variation of its constant 
could resolve some of the standard model problems like G, the constant A is 
a gravity coupling and both should therefore be treated on an equal footing. 

The generalized Einstein's theory of gravitation with time dependent G and A 
has been proposed by Lau [7] . The possibility of variable G and A in Einstein 
theory has also been studied by Dersarkissian [8]. This relation plays an impor- 
tant role in cosmology. Berman [9] and Sistero [10] have considered the Einstein 
field equations with perfect fluid and variable G and for Robertson- Walker line 
element. Kalligas et al. [11] have studied FRW models with variable A and G 
and discussed the possible connection with power-law time dependence of G. 
Abdussattar and Vishwakarama [12] presented R-W models with variable A 
and G by admitting a contracted Ricci collineation along the fluid flow vector. 
Recently some of us and others have studied cosmological models with vari- 
able G and A in a diversified fields [T3lfTlfT5TT6lfTTlfT8^^ Thus 
the implication of time varying A and G are important to study the early 
evolution of the universe. 

In the present paper a five dimensional Kaluza-Klein cosmological model is 
considered with variable G and A which provides an expanding universe. It 
is found that the gravitational constant decreases with time as G ~ £~( 1 ~ a; ) 
whereas the cosmological constant decreases as A ~ t~ 2 , A ~ (R/R) 2 and 
A ~ R/R where u is the equation of state parameter, R is the scale factor. 
However, we have come across with some awkward situations in connection to 
the deceleration parameter and the age of the Universe as far as the scenario 
of present accelerating Universe is concerned. In the present investigation, our 
approach is similar to that of Collins et al. [25] for solving the Einstein field 
equations. The paper is organized as follows: In Sections 2 and 3 we present 
the basic field equations governing the models and their solutions. A critical 
discussion and conclusion are provided in Section 4. 



2 The Einstein field equations for the cosmological model 

We consider five dimensional Kaluze-Klein space-time given by 



ds 2 = dt 2 - A 2 (dx 2 + dy 2 + dz 2 ) - B 2 d^ 2 (1) 

where A and B are function of the temporal coordinate t only. 

The universe is assumed to be filled with distribution of matter represented 
by energy-momentum tensor for a perfect fluid 

Tij = (p + p)ViVi - pg {j (2) 

where p is the energy density of the cosmic matter and p is its pressure, Vi is 
the unit flow vector such that ViV 1 = 1. 

We assume that the matter content obeys an equation of state 

p = up, < u> < 1. (3) 



The field equations are those of Einstein but with time dependent cosmological 
and gravitational constants and are given by Weinberg 



R ij - -Rg ij = 8nGT lj - Ag ij . (4) 

The spatial average scale factor R(t) is given by 



R 4 = A 3 B (5) 

and average volume scale factor V = -R 4 . 

The average Hubble parameter H may be generalized in anisotropic cosmo- 
logical model as 



H = -X = Y = -[ — + — (6) 

4 4 \A BJ K J 

where X = logV and Y = logR and an over head dot denotes ordinary 
differentiation with respect to cosmic time t. We also have 



H=^(H 1 + H 2 + H 3 + H,) (7) 

where H i = H 2 = H 3 = A/A and H4 = B/B are Hubble's factor in the 
directions of x, y, z and ip respectively 



For the line element (1) the Einstein field equations (4) yield the following 
equations 

oA oAd , .. 

+ -T7T = ^Gp + A > ( 8 , 



A 2 AB 



2A B 2AB A 2 

-A+B + ^B+T* = - 8 * Gp + A > (9) 



3 4 3 4 2 

— + -^ = -MG V + K. (10) 

If we use the equivalent energy conservation of general relativity by taking the 
covariant derivative of Einstein field equations (4) we find that 



N [3A B\ G A 

? + (p+ri( T +g)+ g P + ^ = o. (n; 

This provides [T7K2] 



N f3A B\ , . 

p + (p + pA- a - + b) =0 (12) 

as well as 



A = -8nGp (13) 

which will be satisfied by the solutions of the field equations as can be verified 
easily. 



3 The general solutions to the field equations 

The equations (3), (8) - (10) and (12) are five independent equations with six 
unknowns A, B, p, p, G and A. Hence to get a realistic solution we assume 
that the expansion scalar 9 in the model is proportional to the shear a |21j . 
This condition leads to the relation between metric potential [25] 

B = A n (14) 



where n(^ 1) is a constant. 

From equations (9) and (10), we have 

A B ki 



A B A 3 B 



(15) 



where k\ is a constant of integration. Using equation (14) in equation (15) and 
then integrating it we get the lime element (1) as 



ds 2 = dt 2 



, ,n + 3) 

fciTj -{t + k 2 (n + 3) 



-i 2/(n+3) 



1-n 



h^^lt + k 2 (n + 3) 



(dx 2 + dy 2 + dz 2 ) 

2n/(n+3) 



# J 



(1-n) 
By suitable changes of constants the line element (1) can be written as 



(16) 



ds 2 = dt 2 - {at + b) 2/(n+3) (dx 2 + dy 2 + dz 2 ) - {at + b) 2n/(n+3) difj 2 (17) 

where a — k\{n + 3)/(l — n) and b = ^(n + 3) are two positive constants. 

For the model (17), the spatial volume V, matter density p, pressure p, grav- 
itational parameter G, cosmological parameter A are given by 

V = R A = (at + b), (18) 
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3a 2 (n + l) 



4ttA; 2 (1 + u)(n + 3) 2 {at + b)^'^ ' 



(21) 



A 



3(l-w)a 2 (n+l) 
[l + u){n + 3) 2 {at + b) 2 ' 



(22) 



The expansion scalar 9 and shear a are 



e = m= . - (23) 

(at + 6) ' y ' 



(24) 



(at + 6) 
where a is a positive constant. 

The critical density p c , vacuum density p\ and density parameter Q are given 
by 



3^ _ (l+a;)(n + 3) 2 

^ c 87rG 32A; 2 (n + l)(at + 6) (!+-)' l ' 



_ A _ fc 2 (o;-l) 
PA "87rG"2(at + 6)(i+-)' l ° J 



= 87rCp = (n+1) 

3# 2 (l+w)(n + 3) 2 ' l J 



Keeping in mind the form of the solution (18), i.e. R 4 = (at + 6), one can 
express all the above parameters as function of the scale factor R as well 
which we will follow in the following two sub-cases. 



3.1 Special case studies 

3.1.1 The radiation era 

For the radiation epoch, we have u = 1/3. In this case 

p = k 2 R- 4 '\ (2* 



_ 9a 2 (n + 1) 9/ , „ x 

G = \ > R - 2/ 3 , 29 



A = -y^A, (30) 

2(n + 3) 2 ' V ; 



Pu = -^R~ 4/ \ (31) 



(n + 3) 2 4/3 

Pc " W^TT) ■ (32) 



We observe that in this phase matter density is three times the vacuum density. 
The cosmic density parameter Q in this phase is 



a = ^±4 03) 

4(n + 3) 2 V ; 



3. l.£ The dust epoch 

For the phase dominated by dust matters, we have uu — 0. In this epoch 

p = k 2 R~\ (34) 



G= ^ W + j ir 1 , (35) 

47rfc 2 (n + 3) 2 ' v ; 
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8^ = \ < 39 » 



4 Discussion and conclusion 



In the present phenomenological study five dimensional Kaluza-Klein cos- 
mological models with varying G and A are obtained. The models obtained 
present an expansion scalar 9 bearing a constant ratio to the anisotropy in the 
direction of space like unit vector A' and with equation of state p = up, < 
u < 1. If we inspect the solution set then it can be observed from the equation 
(18) that the spatial volume V is zero at t = t where t = —b/a and expan- 
sion scalar 9 is infinite at t — t which shows that the Universe starts evolving 
with zero volume and infinite rate of expansion at t — t . The scale factors also 
vanish at t — t and hence the model has a point type singularity at the initial 
epoch. Since a/9 is constant, as evident from the solutions (23) and (24), the 
anisotropy does not die out asymptotically. The model admits a negative A 
unless u = 1 and these type of models for the Universe with negative A are 
available in the literature [27] • We observe from the general solution set that 
for the stiff fluid (a; = 1) the model reduces to standard Kaluza-Klein type 
Zel'dovich universe [28J with A = 0, G = constant and p = p ~ t~ 2 for the 
limiting case b — > (Figs. 1 & 2). At t = t limit all the parameters G, A, 
a 1 and 9 are infinite in value, whereas all these become zero as t — > oo. The 
ratio between cosmic matter density parameter Q m (= 8irGp/3H 2 ) and cosmic 
vacuum-energy density parameter ^a(= A/3H 2 ) scales as 



^ = -A_ = ^ (40) 

tt A 8nGp 2 ' K J 

whereas cosmic vacuum-energy density parameter can be given as 



J_ g(n + l)(o;-l) 

A " 3H 2 ~ (n + 3) 2 (l + u;)- [ } 



The above two expressions related to cosmic matter and vacuum density pa- 
rameters do vanish for the stiff fluid condition u = 1. 

If we look at the solution expressed in the equation (21) then we can observe 
that in general G ~ t^ 1- ^. Thus for the dust case (to = 0) this case is 
similar to the result of Dirac [I], Canuto et al. [2|3] and Beesham [4]. On 
the other hand, equations (22) and (23) at once imply that A ~ H 2 for any 
constant values of to and n. Again, from the solutions (18) and (22) one can 



find A ~ R/R. On the other hand for the limiting case b — > the equation 
(22) yields the relation A ~ t~ 2 . Thus we obtain A ~ H 2 = (R/R) 2 , A ~ 
R/R and A ~ t~ 2 which is in accordance with the main dynamical laws 
one finds in literature proposed for the decay of A. Note that we find all 
these phenomenological relations for A directly from the solution set which 
generally consider as ad hoc basis in cosmological research. The dynamical 
law A ~ (R/R) 2 , has been proposed by Carvalho et al. [29] and considered by 
Salim and Waga [30], Arbab and Abdel- Rahman [31], Wetterich [32], Arbab 
[33] and very recently by Ray and his collaborators as well as Pradhan and 
his collaborators [T5l[T6 .34,18,35 36]. The decay law of the form A ~ R/R 
has been considered by Arbab [37], Khadekar et al. [H], Ray et al. [15] . Ray 
and Mukhopadhyay [31] and Ray et al. [35]. The dynamical law A ~ f 2 has 
been considered by several authors, e.g. Bertolami [38], Berman and Som |39j, 
Berman [ID] , Beesham [4T] , Pradhan et al. [T3] to mention a few. We observe 
that |A| decays faster than G where as p, p, p v , p c scales R~^ +1 ^ as can be 
seen from the Figs. 1 - 4. In the model we see that the quantity G satisfies the 
condition for a Machian cosmological solution i.e. Gp ~ H 2 [42j. 



However, as far as the scenario of present accelerating Universe is concern 
our models seems suffer a lot in connection to the different physical param- 
eters, such as the deceleration parameter and the age of the Universe. The 
deceleration parameter q for the present model is q = 3 and hence does not 
correspond to the observational result go ~ —0.5 [4"3"f4"4"] . However, we have 
observed that for the scale factor of the form a(t) = t n we can have a negative 
deceleration parameter when n > 1. In our present case n = 1/4 which yields 
q = +ve. Here the age of the universe is given by t = \Hq 1 which is also 
differs from the present estimate to = Hq 1 « 14 Gyr |4"5|46|I47||4'8] . This value, 
therefore, obviously suffers from the low age problem [19] as even the oldest 
globular clusters show a value which is t ~ 11 Gyr [50f5T] . In this connection 
it is argued that dark energy models without A suffer from low age problem 
[52f4"9] though even in the presence of A we are not in a hopeful position. One 
of the reasons of these discrepancies may be due to the constant nature of the 
equation of state parameter u which can not accommodate with the present 
accelerating phase of the Universe. However, this needs further investigation. 
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Fig. 1. Variation of matter density with time for the given values of u under the 
constraint k2 = 0.8, a = 1 and b = —2. 
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Fig. 2. Variation of pressure with time for the given values of u>. 
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Fig. 3. Variation of time-dependent gravitational constant with time for the given 
values of u. 
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Fig. 4. Variation of time-dependent cosmological constant with time for the given 
values of to. 
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